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Resumen 


La cuenca del río Dróme, localizada en el sureste de Francia, ha experimentado en 
las décadas recientes periodos recurrentes de baja disponibilidad hídrica y crisis 
severas de escasez de agua en el sector agrícola. El índice normalizado de lluvia 
(SPL, por sus siglas en inglés, Standardized Precipitation Index) se ha utilizado 
para caracterizar la sequía en términos de su frecuencia, intensidad y duración. En 
esta investigación se propone la aplicación de la distribución gamma generalizada 
extendida para calcular los índices SPI e identificar la sequía estacional en la cuenca 
del río Dróme con base en el análisis de registros históricos de lluvias de la estación 
meteorológica de Montelimar, Francia. Los parámetros de la distribución se 
estimaron numéricamente por el método de máxima verosimilitud. La bondad de 
ajuste del modelo propuesto fue probada con el estadístico de Kolmogorov-Smirnov 
usando bootstrap paramétrico. Los índices SPI fueron calculados para un año 
agrícola y para las épocas de riego de primavera-verano y verano. La distribución 
gamma generalizada extendida presentó un buen ajuste a las frecuencias de lluvias 
acumuladas observadas. La evolución de los índices SPI permite constatar una 
tendencia negativa de la lluvia estacional en primavera y verano en las dos últimas 
décadas en la región bajo estudio. 


Palabras clave: sequía meteorológica, variabilidad climática, máxima verosimilitud, 
bootstrap paramétrico. 


Introducción valle, los cultivos más importantes son maíz, 
girasol, trigo, hortalizas y plantas aromáticas. 
En los años recientes, el déficit hídrico aumentó 
principalmente debido a la variabilidad de la 


lluvia y al incremento de la superficie de riego 


En las últimas dos décadas, la ocurrencia de 
años secos ha sido cada vez más frecuente en la 
cuenca del río Dróme, localizada en el sureste 


de Francia. En esta región mediterránea, la 
irrigación es un factor esencial para garantizar 
la producción agrícola y desde 1997 cuenta 
con un plan de ordenamiento y gestión del 
agua denominado SAGE (por sus siglas en 
francés, Schéma d'Aménagement et de Gestion 
des Eaux), para disminuir el impacto negativo 
sobre el río Dróme, inducido por un incremento 
de la demanda y de la contaminación. En este 


(Gleyses y Rieu, 2004). El déficit hídrico para 
uso agrícola en la cuenca es de 2 a 9 millones 
de m?; un año sobre cuatro; los caudales del río 
presentan un coeficiente de variación (CV) de 
62% en julio y de 96% en septiembre (González- 
Camacho ef al., 2008). 

Con el propósito de caracterizar las ten- 
dencias y los periodos de sequía se calcularon 
los índices normalizados de lluvia (SPI, por sus 
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siglas en inglés, Standardized Precipitation 
Index) propuesto por McKee et al. (1993). 
Estos índices han sido utilizados por su 
simplicidad de interpretación para el moni- 
toreo de sequías a diferentes escalas de tiempo 
(McKee et al., 1995). El cálculo de los SPI se 
realiza comúnmente con base en el análisis de 
frecuencias de lluvias mensuales acumuladas 
para diferentes intervalos de tiempo y en el 
uso de la distribución gamma de dos pará- 
metros, los cuales se estiman mediante una 
aproximación analítica de la función de má- 
xima verosimilitud; luego las probabilidades 
teóricas de lluvia acumulada son normaliza- 
das mediante la función de cuantiles de la 
distribución normal estándar, que proporciona 
valores relativos entre -3 y +3 (Thom, 1958). 
Mediante esta transformación, los valores 
positivos de los SPI indican cantidades de llu- 
via superiores a la mediana y valores negativos 
corresponden a cantidades de lluvia inferiores 
ala mediana. Guttman (1999) propone el uso de 
la distribución Pearson tipo lll o gamma de tres 
parámetros como una función más versátil. 
Vicente-Serrano (2006) utiliza esta distribución 
para realizar un análisis espacio-temporal 
de las sequías en la península ibérica. Sin 
embargo, la distribución gamma generalizada 
extendida de tres parámetros descrita por Lee 
y Wang (2003) es una función aún más general 
y más versátil para ajustar los datos de lluvia 
acumulada. Dos dificultades para utilizar 
modelos más generales y más flexibles son los 
procedimientos de estimación de sus paráme- 
tros y de la prueba estadística de bondad de 
ajuste de los mismos, que en general requieren 
soluciones numéricas computacionales. Sin 
embargo, la disponibilidad de tecnología 
informática de gran capacidad de cálculo y 
bajo costo hace factible el uso de modelos más 
complejos para resolver problemas prácticos, 
como lo es el cálculo de los índices SPL. 

Los índices se calculan para diferentes 
intervalos de tiempo. Los SPI-6 y SPI-12 
permiten caracterizar el comportamiento 
semestral y anual de la lluvia a mediano plazo; 
el índice SPI-3 proporciona una comparación 


de la lluvia para intervalos trimestrales de la 
serie histórica. Esto es, un SPI-3 agosto compara 
la lluvia total de junio, julio y agosto de un 
año específico con las lluvias totales de junio 
a agosto de los demás años. El índice SPI-3 
refleja las condiciones de humedad a corto 
plazo y proporciona una estimación estacional 
de la lluvia. En México, los índices SPI han sido 
utilizados para caracterizar algunas regiones 
áridas y semiáridas del norte de México (entre 
otros, Velasco et al., 2004). 

Por lo anterior, en la presente investigación 
se propone aplicar la distribución gamma ge- 
neralizada extendida para calcular los índices 
SPI y caracterizar la sequía estacional en la 
cuenca del río Dróme con base en el análisis 
de registros históricos de lluvias de la estación 
meteorológica de Montélimar, Francia. Los 
parámetros de la distribución se estiman 
numéricamente por el método de máxima 
verosimilitud y se utiliza una técnica bootstrap 
paramétrica para calcular los valores críticos 
del estadístico de la prueba de Kolmogorov- 
Smirnov. 


Metodología 
Registros históricos de lluvia 


Las series de lluvias mensuales se obtuvieron 
a partir de registros diarios de lluvia (1960- 
2006) de la estación de Montélimar (latitud 
44.35” N, longitud 4.44” E, altitud 73 msnm), 
la cual es representativa de la zona de estudio. 
Se generaron las series de lluvia acumulada 
para los siguientes intervalos de tiempo: 
anual (octubre-septiembre), que corresponde 
al año agrícola; semestral (marzo-agosto), que 
corresponde a la época de riegos de primavera- 
verano, y trimestral (junio-agosto) para la 
época de riegos de verano (periodo estival). Se 
calcularon los estadísticos descriptivos media, 
mediana, desviación estándar y el coeficiente 
de variación (CV = desviación estándar / media) 
expresado en por ciento, para caracterizar las 
series de lluvias analizadas. 
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Distribución gamma generalizada extendida 


Se asume que la distribución de frecuencias 


empíricas de lluvias para los diferentes 
intervalos de tiempo estudiados puede aproxi- 
marse con la función de densidad gamma 
generalizada extendida (GGE) descrita como 
sigue (Lee y Wang, 2003). 

Una variable aleatoria X se dice que tiene 


distribución GGE (a, A y) si su función de 


densidad está dada por: 
lafy ay -1 o 
fío Ay) =.—8 exp l-y(x) L, 
Pm) 
x>0, y>0, 2>0,100€R (1) 


donde a. y y son parámetros de forma; A es 
un parámetro de escala, y (y) es la función 
gamma, dada por T(y) = ' epi-tHat 
(Abramowitz y Stegun, 1965). 

La distribución GGE es bastante flexible, 
tiene como casos especiales la distribución 
gamma (a = 1), la distribución exponencial (0. 
= y = 1), Weibull (a. > 0; y = 1) y la lognormal 
cuando y > %, algunas de las cuales son 
ampliamente utilizadas en hidrología. 

La función acumulativa de probabilidades 
de X está dada por (Lee y Wang, 2003): 


Elx;0,A y) =P(X< x) 


ls Tu(y0x)7), a.<0 
= : 0) 
Ty (10) 7 1) a>0 


donde Pu(x,y)= 0 [*TopEgat es la fun- 
Y 


ción gamma incompleta (Didonato y Morris, 
1986). 

Los parámetros al, A y y se estiman mediante 
la maximización numérica de la función 
de verosimilitud. Suponiendo que se tiene 
la muestra X,,...,X, de la distribución GGE 


(a,2,y), la función de verosimilitud que debe 
maximizarse es: 


n 


My. .., L)= TTsa(x;: 0 A,y) 


¡i=l 


(0, Y (3) 


O bien una función monótona de la misma: 


1(0, yx, a a) =l0gL(0,»y|x. ' $e) 


= 


i 


log f., (x;;0, Ay) 
1 


La maximización de la ecuación (3) se 
realiza mediante técnicas numéricas, cuyos 
algoritmos son ilustrados por Stacy y Mihram 
(1965); igualmente, la maximización puede 
realizarse con el procedimiento LIFEREG del 
paquete estadístico SASO (SAS, 2008). 


Prueba de bondad de ajuste de la distribución 
GGE 


Dado x,,...x, (i.e. los datos de lluvia acumula- 
da), se desea saber si el modelo GGE aproxima 
adecuadamente la distribución de frecuencias 
empíricas de lluvia. La prueba de bondad 
de ajuste consiste en probar las hipótesis 
siguientes: 


Bo :EJx)=F,(x;0,A,y) vs 
(4) 
Hi: E,(x)+ E x;0, A y) 


Existen varios estadísticos de prueba para 
probar el juego de hipótesis en la ecuación (4), 
siendo uno de los más conocidos y usados el 
de Kolmogorov-Smirnov (Massey, 1951). Dado 
z, =Flx,), x, < ... <x,, donde F(e) es la función 
acumulativa de probabilidades de GGE (los 
parámetros desconocidos son reemplazados 
por sus estimadores de máxima verosimilitud), 
el estadístico de prueba de Kolmogorov-Smir- 
nov es: 


Dar =máx(D*, D”) (5) 
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donde: 


i ¡IN 
+ A E 
D a E máx cisn[3= me ] 


La regla de decisión es rechazar H, si D.,,. 
> D, (a) con a e (0,1) el nivel de significan- 
cia de la prueba. Los valores críticos de la 
prueba de Kolmogorov-Smirnov son válidos, 
siempre y cuando los parámetros de la 
distribución F sean completamente conocidos; 
en este estudio, como éstos se estiman por 
máxima verosimilitud, ya no son aplicables, 
por lo que se propone utilizar una versión 
bootstrap de la misma para obtener los valores 
críticos mediante el procedimiento bootstrap 
paramétrico siguiente: 


1. Dadas las observaciones x.,..., x,, Calcular 
los estimadores de máxima verosimilitud 


de 0 = (0,2,y), esto es, Ó= (ó, A 1) 


eo Utilizando la ecuación (5) y 

los estimadores 6 = do, Y a 

Il. Generar una muestra bootstrap de 
tamaño n de F,= (a,4,7) el algoritmo 
de generación de números aleatorios es 
descrito por Stacy y Mihram (1965). 

Il. Dada la muestra bootstrap generada 
previamente, calcular los estimadores 
de máxima verosimilitud de O = (0, A,y), 
¡.e. Q= (5,4, : 

III. Calcular D', usando los estimadores 
9= (ó, de 7) y la muestra bootstrap. 


3. Repetir los pasos l), II) y III) mil veces para 


obtener Do; ¡;j=1 ...11000. 


4. Fijar el nivel de significancia de la prueba, 


por ejemplo a.=0.05 y obtener el valor de la 
constante crítica D,, ,.. (0.05) como ay 
donde ct ;j=1, ... 1 000 denota los 


valores ordenados de D* 


calc* 


5. Se rechaza H, al nivel de significancia Ol si 


D >D' 


calce calc,(950)* 


El bootstrap paramétrico se implementó 
en un programa IMLQ de SASO (SAS, 
2008), que puede obtenerse solicitándolo a 
cualquiera de los autores o en la dirección: 
http: / /174.123.176.26/SAS_macro/,  recupe- 
rando el archivo gamma_generalizada.sas. 


Cálculo de índices normalizados de lluvia SPI 


Dado que la distribución GGE no está definida 
parax=0y que una distribución de lluvia puede 
tener valores nulos, se utiliza una distribución 
GGE ponderada, que se expresa por: 


H(x)=3+(1-3)Fx(x;0,2,y) (6) 


Donde q es la probabilidad de tener un 
valor igual a cero. La probabilidad q puede ser 
aproximada según Thom (1958) por el cociente 
q = m/n, donde m es el número de datos con 
lluvia acumulada igual a cero y n es el número 
total de valores acumulados de lluvia. Luego 
H(x) se transforma, por medio de la función de 
cuantiles de la distribución normal estándar a 
una variable Z reducida, denominada índice 
SPI; una aproximación de esta función está 
dada por Abramowitz y Stegun (1965): 


2 
otaHtarI. e 0<H(x)<0.5 
1+d, t+ d, t*+d,t 
Z=SPlI= 
2 
Cote tb 0.5<H(x)<1.0 
l+ d, t+ d,tó+d,t 
(7) 
donde: 
1 , 
In Si si0<H(x)<05 
H(x) 


si0.5<H(x)<1.0 


de 
l (1-H()) 


Cy =2.515517; c, =0.802853; c7 = 0.010328; 
d¡ =1.432788;d, =0.189269; dz =0.001308; 
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Cuadro 1. Clasificación de los índices normalizados de 
lluvia SPI (McKee ef al., 1993). 


SPI Descripción 

> +20 Muy húmedo 
de+l9Na 15 Húmedo 
de +1.49 a+1.0 Moderadamente húmedo 

+0.99 a -0.99 Normal 

de -1.0 a -1.49 Moderadamente seco 
de -1.5 a -1.99 Seco 

<-2.0 Muy seco 


La clasificación de los SPI calculados me- 
diante la ecuación (7) para caracterizar las 
lluvias acumuladas se presentan en el cuadro 
1. 


Resultados y discusión 


Distribución de la lluvia acumulada de 
octubre a septiembre (1960-2006) 


La distribución interanual de la lluvia acumu- 
lada de octubre a septiembre en el periodo 
1960-2006 tiene una variabilidad interanual 


(expresada por el CV) de 18%, con una media 
de 919+/-168 mm y una mediana de 907 mm 
(figura 1). El año más lluvioso se presentó en 
el periodo de octubre de 1976 a septiembre de 
1977 con 1 388 mm, y el más seco en el periodo 
de octubre de 1989 a septiembre de 1990 con 
544 mm. Entre 1986 y 1990 se presentaron, en 
cuatro años de cinco, lluvias inferiores a 800 
mm, provocando periodos de déficits hídricos 
en la cuenca del río Dróme. 


Distribución de la lluvia acumulada en 
primavera-verano y en periodo estival 


La distribución interanual de la lluvia en el 
periodo de primavera-verano (de marzo a 
agosto) muestra una fuerte variabilidad (CV 
= 30%). El promedio de lluvia acumulada es 
de 387+/-116 mm y la mediana es de 367 mm 
(figura 2). La lluvia acumulada durante este 
periodo representa aproximadamente 40% de 
la lluvia anual de octubre a septiembre. 

La distribución interanual de la lluvia en 
el periodo estival (de junio a agosto) tiene 
un CV = 47%. La lluvia media estival es de 
168+/-82 mm y la mediana es de 146 mm; 


Media = 920 mm CV=18% 


Lluvia (mm) 


bl . E 7 ER | Al ñ GS 
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Año 


Figura 1. Distribución de la lluvia anual acumulada de octubre a septiembre (1960-2006) en Montélimar, Francia. 
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Figura 2. Distribución de la lluvia acumulada de marzo a agosto (1960-2006) en Montélimar, Francia. 
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Figura 3. Distribución de la lluvia acumulada de junio a agosto (1960-2006) en Montélimar, Francia. 
el valor mínimo de 29 mm sucede en 1989 y Estimación de los parámetros de la función de 
el máximo de 422 mm ocurre en 1992 (figura densidad GGE 
3). En este periodo, el riego es esencial para 
cubrir las necesidades hídricas de los cultivos. Los parámetros de la función GGE estimados 
Asimismo, de 1990 a 2006 se observan pe- por máxima verosimilitud para las diferentes 
riodos estivales menos lluviosos que en el series de lluvia estudiadas se presentan en el 
periodo 1960-1989. cuadro 2. 
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Cuadro 2. Parámetros de máxima verosimilitud de la 
distribución gamma generalizada extendida, de lluvias 
acumuladas de doce, seis y tres meses. 


Parámetros* de la 
distribución GGE 


Lluvia SPI-12 


z 0.28173 | 0.0010998 | 389.923 
septiembre 
Lluvia SPI-6 agosto | -0.61714 | 0.0027651 | 32.322 
Lluvia SPI-3 agosto 1.27854 | 0.0057849 2.624 


* a y yson parámetros de forma, y A es un parámetro de 
escala de la distribución gamma generalizada extendida 
(GGE) (Lee y Wang, 2003). 


Prueba de bondad de ajuste de la función de 
distribución GGE 


A partir de la función de densidad GGE 
(ecuación (1)) se calcularon las distribuciones 
de probabilidades teórica y empírica. La 
prueba de bondad de ajuste fue realizada 
utilizando la versión bootstrap de la prueba 
de Kolmogorov-Smirnov para la distribución 
GGE. Para un nivel de probabilidad a = 0.05, 
las distribuciones teóricas de probabilidades 
de lluvia acumulada no son diferentes de la 
distribuciones de frecuencias de lluvias para 
las series estudiadas. Los resultados de las 
pruebas estadísticas para los tres intervalos de 
tiempo estudiados se presentan en el cuadro 3. 

En la figura 4 se muestra la distribución 
del estadístico Dmáx de Kolmogorov-Smir- 
nov obtenido con el procedimiento bootstrap 
paramétrico, para los datos de lluvias acu- 
muladas de octubre a septiembre (1960-2006). 


0.11 0.13 0.15 0.17 0.19 


Dmáx 


Figura 4. Distribución bootstrap de Dmáx SPI-12 septiembre, 
Montélimar, Francia. 


Además, en la figura 5 se muestra la bondad 
de ajuste entre la distribución empírica 
de frecuencias de lluvia acumulada y la 
distribución GGE ajustada para la serie de 
lluvias acumuladas de octubre a septiembre 


(1960-2006). 
Índices normalizados de lluvia SPI 


La distribución interanual de los índices SPI-12 
septiembre (figura 6) muestra una disminución 
de años muy húmedos a húmedos (SPI entre 
1 y 2.5) en el periodo 1990-2006 con respecto 
al periodo 1960-1989 (SPI entre O y 1). En 
cambio, para los mismos periodos, los años son 
menos secos (SPI entre O y -1.0) en el último 
periodo que en el anterior (SPI entre -0.5 y 
-1.5). El año 1990 fue extremadamente seco 
(SPI = -2.5). Igualmente, se observa en todo el 
periodo estudiado una frecuencia alta de años 
normales (SPI entre -0.99 y +0.99). Sin embargo, 


Cuadro 3. Prueba de bondad de ajuste de Kolmogorov-Smirnov para la distribución gamma generalizada extendida. 


Lluvia SPI-12 octubre-septiembre 


Lluvia SPI-6 marzo-agosto Lluvia SPI-3 junio-agosto 


» 0.0513 0.1047 0.1065 
Dosítico DOOÍSHAp 0.19 0.21 0.24 
H,; GGE = empírica No rechazo No rechazo No rechazo 


Se rechaza H,: GGE = empírica; si D,,,, > D, 


crítico 


bootstrap con un nivel de probabilidad a = 0.05. 
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Figura 5. Bondad de ajuste de la distribución GGE a la distribución empírica de frecuencias de lluvias, SPI-12 septiembre, 


en Montélimar, Francia. 
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Figura 6. Distribución interanual de los SPI-12 septiembre (1960-2006) en Montélimar, Francia. 
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Figura 7. Distribución de los índices SPI-6 agosto (1960-2006) en Montélimar, Francia. 


la variabilidad climática es más fuerte en el 
periodo anterior a 1990 (SPI entre -1.5 y 2.0) 
que después de 1990 (SPI entre -0.5 y 1.0). 

La distribución interanual de los índices 
SPI6 agosto (figura 7), que toma en cuenta 
las lluvias de marzo a agosto, describe una 
frecuencia de años más secos mayor después 
de 1990 (SPI entre -1.0 y -2.0) respecto al 
periodo 1960-1989. De 1990 a 2006 ocurrie- 
ron siete años secos (SPI entre -1.0 y 2.0), 
mientras que de 1960 a 1989 ocurrieron tres 
años secos (SPI entre -1.0 y -2.0). En cuanto 
a los años húmedos, la frecuencia de años 
húmedos es menor después de 1990 (SPI entre 
+1.0 y +2.5) que antes de 1990. Igualmente 
puede observarse que estacionalmente, 2003, 
2004 y 2005 se consideran años moderada- 
mente secos y secos, mientras que para un 
periodo anual estos años fueron normales. 
Esto permite explicar que en los últimos 
años los riegos en primavera y verano sean 
cada vez más frecuentes. Asimismo, para los 
cultivos de primavera se requiere un riego de 
pre-siembra para asegurar la germinación de 


los cultivos de verano (González-Camacho et 
al., 2008). 

La distribución interanual de los índices 
SPL3 agosto (figura 8), que toma en cuenta las 
lluvias acumuladas de junio a agosto, muestra 
una mayor variabilidad de los índices de 
lluvia (SPI entre -2.5 y +2.5). El año 1989 fue 
seco en el verano (SPI=-2.5). Globalmente, en 
el periodo de riegos de verano en los últimos 
veinte años, 14 años sobre 20, las lluvias 
acumuladas de junio a agosto estuvieron por 
abajo de la media, de los cuales siete años se 
consideran como secos o muy secos. Estos 
índices muestran que las lluvias estacionales 
presentan una tendencia negativa en los 
últimos 15 años. Desde 1989 se observa que 
la ocurrencia de años secos es cada vez más 
frecuente. Para la zona estudiada, 1989 fue el 
más crítico, ya que los tres índices estudiados 
clasifican el año de 1989 como seco o muy 
seco. 

El análisis de los índices SPI permite 
constatar un aumento de las frecuencias de 
periodos secos o muy secos en los últimos años 
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Figura 8. Distribución de los índices SPI-3 agosto (1960-2006) en Montélimar, Francia. 


en la cuenca del río Dróme. De 1960 a 1989, la 
presencia de años moderadamente húmedos 
es más frecuente (SPI entre +1.0 y +1.5); mien- 
tras que en los últimos años, la frecuencia 
de años moderadamente secos a muy secos 
son más frecuentes (SPI entre -1.0 y -2.5). 
Asimismo se observa una tendencia negativa 
más acentuada de las lluvias estacionales 
SPI-6 agosto y SPI-3 agosto (figuras 7 y 8) 
en relación con las lluvias anuales SPI-12 
septiembre (figura 6). 


Conclusiones 


El cálculo de los índices normalizados de 
lluvia SPI mediante la distribución gamma 
generalizada extendida y el análisis de la 
distribución de la lluvia estacional en la cuenca 
del río Dróme, Francia, en periodos, donde 
la práctica de la irrigación es esencial para 
garantizar la producción agrícola de esta región 
permiten obtener las siguientes conclusiones: 


La 
extendida, cuyos parámetros se estimaron 


distribución gamma generalizada 
por máxima verosimilitud, presentó un 
buen ajuste a la distribución empírica de 
frecuencias de lluvias observadas para 
los intervalos de tiempo estudiados (año 
agrícola, primavera-verano y verano), en la 
cuenca del río Dróme. 

El procedimiento bootstrap ¡paramétrico 
propuesto permitió calcular los valores 
críticos de la prueba de bondad de ajuste de 
Kolmogorov-Smirnov para los diferentes 
de El 
procedimiento de estimación y prueba de 


intervalos tiempo estudiados. 
bondad de ajuste de la distribución gamma 
generalizada extendida está disponible en 
la dirección http://174.123.176.26/SAS_ 
macro/, recuperando el archivo gamma_ 
generalizada.sas. 

Los índices SPI-12 septiembre varían entre 
-1 y 1 en los años recientes, por lo cual los 


años se consideran como normales. Sin 
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embargo, estacionalmente, en primavera- 
verano y verano, los índices SPI-6 agosto y 
SPI-3 agosto con valores entre -1.0 y 2.5 
son más frecuentes de 1990 a 2006, lo que 
muestra un incremento en la ocurrencia de 
periodos cada vez más secos en primavera- 
verano y verano en la cuenca del río 
Dróme. 
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Abstract 


GONZÁLEZ-CAMACHO, J.M., PÉREZ-RODRÍGUEZ, P. € RUELLE, P. Estimation of 
standardized precipitation indexes using the extended generalized gamma distribution. Water 
Technology and Sciences (in Spanish). Vol. II, No. 4, October-December, 2011, pp. 65-76. 


The Dróme river basin is located in southern France, where in recent decades the agricultural 
sector has undergone recurrent periods of low water availability and severe water scarcity. 
The Standardized Precipitation Index (SPI) has been used to characterize drought in terms of 
frequency, intensity and duration. In this research, the use of the extended generalized gamma 
distribution is proposed to compute the SPI and to identify seasonal drought in the Dróme 
river basin based on the analysis of historical rainfall data from the meteorological station 
in Montélimar, France. The parameters were numerically estimated using the maximum 
likelihood method. The goodness-offit test of the model was based on the Kolmogorov-Smirnov 
statistic, computed with a parametric bootstrap technique. The SPls were computed yearly 
and for the spring-summer and summer irrigation seasons. The extended generalized gamma 
distribution presented a good fit to observed rainfall frequencies. The development of SPI 
indexes enabled verifying a negative trend for seasonal rainfall in spring and summer over the 
last two decades in the region under study. 


Keywords: meteorological drought, climatic variability, maximum likelihood estimation, 
parametric bootstrap. 
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